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SUMMARY 

In this paper we consider the time slot assignment for a heterogeneous environment in which circuit- 
switched traffic and packet-switched traffic share the same satellite. In constructing a single time 
division multiple access (TDMA) frame for both traffic types, their different characteristics must be 
taken into account. This problem is known to be NP-complete and a couple of heuristic partial 
optimization algorithms have been developed. In this paper, we first provide a theoretical result to 
improve the existing partial optimization algorithms; then a fully optimizing heuristic algorithm is 
presented. Simulation results show that our algorithm provides a much better solution quality than 
existing ones. © 1997 John Wiley & Sons, Ltd. 
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1. INTRODUCTION 

As the demands of satellite cornrnuni cation rapidly 
erow. the natural resources thev use. the RF snec- 
trum and the geosynchronous orbit, are becoming 
highly crowded. The utilization of the above natural 
resources is optimized by employing multibeam 
antennas and satellite-switched time-division mul- 
tiple-access (SS/TDMA) techniques. In an 
SS/TDMA system, a satellite has a number of spot 
beam antennas, each providing coverage for a lim- 
ited geographical zone. The solid-state RF switch 
on-board the satellite provides connections between 
the various uplink and downlink beams according 
to the TDMA frame. The frame is divided into time 
slots. Each time slot represents a particular switching 
matrix configuration, which transmits a certain num- 
ber of packets between the connected uplink and 
downlink beams without conflict. The objective is 
to permit the transmission of a given traffic load, 
and to accomplish this with maximum transponder 
utilization. 

Most previous studies on SS/TDMA time slot 
assignment 1 " 7 * 8 - 9 have focused on a homogeneous 
traffic environment, consisting of either slowly vary- 
ing circuit-switched traffic or rapidly varying packet- 
switched traffic. Recently, Bonuccelli, et al. w exam- 
ined the time slot assignment for a heterogeneous 
environment in which circuit- and packet-switched 
traffic share the same satellite. Clearly, circuit- and 
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packet-switched traffic have vastly different charac- 
teristics and requirements. Circuit-switched traffic 
(static traffic) is usually composed of a continuous 
bit-stream of a predetermined bit-rate and requires 
real-time delivery. Packet-switched traffic (dynamic 
traffic) is bursty in nature and can be delayed. In 
constructing a single TDMA frame for both traffic 
types these differences must be taken into account. 
Bonuccelli, et al. xo proved that this problem is 
NP-complete. 

In the system studied, the switchpoint closures in 
each frame slot are determined in response to the 
changing demand for each frame. A low-capacity 
order- wire facility (see Reference 5) is assumed for 
sending connection requests to the satellite. In one 
possible procedure, each earth station within a zone 
transmits to the satellite the number of requests for 
connection to each zone for circuits and packets. 
For each zone, the onboard scheduler calculates the 
circuit demand for each path and forms a circuit- 
switched traffic matrix. It also calculates the new 
and previously queued packet demand for each path, 
and forms a packet-switched traffic matrix. For each 
frame the scheduler determines a switch assignment 
for each frame slot, first for the circuits, and then 
for the packets. It then allocates particular stations 
to the assignments and broadcasts the schedule to 
the ground. 

Bonuccelli, et al. w suggested a partial optimiz- 
ation approach for this problem. The algorithm in 
Reference 10 starts with a given TDMA frame for 
the transmission of the static circuit-switched traffic. 
As switching conflicts usually prevent full utiliz- 
ation, it is likely that there will be some unused 
slots in the frame. The algorithm attempts to fill 



CCC 0737-2884/97/050197-7 $17.50 
© 1997 John Wiley & Sons, Ltd. 



Received March J 995 
Revised September J 995 



c C 

« -M KIM ANI1 5; KIM V 



198 S.-H. KIM AND S. KIM 



these empty slots with packets. The object is to 
maximize the utilization of these empty slots. This 
was an excellent idea, but the simulation result was 
not sufficient because the slots did not optimize 
fully. 

In this paper, we will show that it is possible to 
generate a TDMA frame for the static circuit- 
switched traffic more widely, so that more packet- 
switched traffic can be scheduled. Based on this 
observation, we can immediately improve the algor- 
ithm studied in Reference 10. We suggest a new 
fully optimizing algorithm. Unlike the algorithm by 
Bonuccelli, et a/., 10 our algorithm generates switch- 
ing matrices for the static traffic and the dynamic 
traffic interactively. Simulation results show that our 
algorithm generates much better solutions than that 
of Reference 10. 

This paper is organized in the following way. The 
problem formulation is given in Section 2. Section 3 
discusses the partial optimization approach. 
Section 4 presents our fully optimizing algorithm. 
Section 5 reports and compares the extensive simul- 
ation test results of our algorithm with those of the 
previous algorithm and Section 6 concludes. 

2. NOTATIONS AND PROBLEM 
FORMULATION 

The total number of zones in the system is N and 
we shall label the zones from 1 to N. A traffic 
matrix R is an NxN matrix with non-negative 
integer entries. Entry r 0 in R represents the traffic 
requirement to be transmitted from zone / to zone 
/ A line of a matrix is one column or row of that 
matrix, and the line sum is the sum of all entries 
on one line. The total amount of traffic in R, £,/(,■, 
is denoted by \R\. 

A switching matrix is a traffic matrix with one 
non-zero entry at most in any line. Thus, a switching 
matrix represents traffic that can be transmitted with- 
out conflict. The length or duration of a switching 
matrix Q, L(Q) is the magnitude of the largest entry 
in Q. The transmission of Q takes L(Q) consecutive 
time slots. A line of switching matrix Q is exposed 
if all entries in that line are zero. A line with one 
non-zero entry is called a covered line. 

A transmission schedule (or time slot assignment) 
of a traffic matrix R is a decomposition of R into 
several NxN switching matrices Q\,Q 2 , • -,£?«, such 
that R = Q x + Q 2 +. . .+ Q m . The total length or trans- 
mission time of a schedule is the sum of the lengths 
of its switching matrices: L = A schedule 

for a traffic matrix R is optimal if it achieves 
maximum transponder utilization among all possible 
schedules for R. It has been shown in References 3 
and 4 that this is equivalent to minimize the sched- 
ule length, L. 

The time slot assignment problem studied in this 
paper can be stated as follows. Given NxN circuit- 
switched traffic matrix S and packet-switched traffic 
matrix D, construct schedules S u S 2 ,...S m and 



D ly D 2 ,' yD my such that the amount of scheduled 
dynamic traffic (Sj^D*! is maximized while the 
following equations are satisfied: 



S k + D k is a switching matrix , 
for k — l,...,m , 

L(S k + D k ) = L(S k \ for k = 1,. . .,m , 
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where T max is the largest line sum of S, and m is 
the number of switching matrices. If T msK is smaller 
than the capacity of one frame, the capacity which 
is not designated for circuit-switched traffic is 
reserved for packet-switched traffic only. 

We shall refer to matrix D as the dynamic traffic 
matrix, to D k as the Ath dynamic switching matrix, 

iu o <u liiu Suxliv u a j. ill/ UlauiA aim iQ O k ao uic nxn 

static switching matrix. Equations (1) and (2) ensure 
that the dynamic traffic is scheduled without any 
change to the static traffic. Equations (3) and (4) 
imply that the static traffic is scheduled optimally. 
Equation (5) ensures that the total amount of 
dynamic traffic scheduled does not exceed the 
traffic requirement. 

We first derive an upper bound, UB, on the 
maximum amount of dynamic traffic that can be 
scheduled. Let /?, and Cj be the line sums of row i 
and column j 9 respectively, of the static traffic matrix 
£. Then T max - R, and T max - C, are the total number 
of empty slots in the ;th row and yth column, 
respectively, of the static switching matrices. Simi- 
larly, let DRi and DCj be the line sums of row i 
and column y, respectively, of dynamic traffic matrix 
D. The following UB then becomes an upper bound 
on the maximum amount of dynamic traffic that can 
be scheduled: 



UB = min- 



2 min ( ^max - RhDRi) > 



N 



2 min(7 max -CjJDCj) 



Simulation tests show that this upper bound is better 
than that in Reference 10. 



3. PARTIAL OPTIMIZATION APPROACH 

Bonuccelli, et al. x0 suggested partial optimization 
algorithms for this problem which optimize the 
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scheduling for D only. These algorithms start with 
a given TDMA frame for the static traffic. Then 
they assign as many packets as possible to empty 
slots of the frame. Among their algorithms, the 
local-optimal algorithm is more practical in the sense 
that it has relatively low complexity and provides 
relatively good solutions. 

The partial optimization approach provides a sol- 
ution quickly. However, its solution quality is not 
guaranteed. We will provide a theorem which shows 
a way of improving the quality of the partial optimal 
solution. To prove it, we start with a lemma. 

Lemma 1. S and D are switching matrices such 
that S + D is a switching matrix and L(S + D) = L(S). 
Suppose that S is decomposed to {S 1 ,...,^}, such 
that S=2,US' <mdL(S) = !% l L(S l y Then there is a 
decomposition of D {£>\... such that 

D - SJLi^, & + is a switching matrix and 
LOS") = L(S* + D 1 ) for /= 1,...,?. 

Proof. Let us consider the case where q — 2. Sup- 
pose that {D x jy*} is a decomposition of D, such 
that D = D' +D 2 and D 1 and D 2 are switching mat- 
rices. Since S + D is a switching matrix, S 1 + D 1 is 
also a switching matrix for / = 1,2. Hence 

KS 1 + &) = m3x{L(S , )M&)\ for '=1,2. Now let 
us construct a decomposition as follows: dj, 

= mm{d ij9 L(S 1 )} and 



4 = d,-4 = 



m 



m 



{ 



^-L(S J ), if d,>L(S l ), 
0 , otherwise , 



where d 0 ,d^j are the (iV)th elements of matrices D 
and ID 1 , respectively. Then L(D X ) < L(S ] ). And 
LCD 2 ) is given by 



L(D) -L(S l ) 9 if L(D) > L(S l ) , 
0 , otherwise . 



Since L(D) < L(S), L(Z) 2 ) < - L(S X ) = L(S 2 ). 
Hence + D 7 ) = max ),£(£>')} = for 
/=1,2. This completes the proof for q = 2. The 
proof for q > 2 immediately follows by successively 
applying the result for = 2. 

Theorem 1. Let {5 , l ,...,S m } be given switching mat- 
rices for the static traffic. Suppose that each S k is 
decomposed to switching matrices {S k , - ,S% k } such 
that 

(a) S k = l% S[., 

(b) L(S k ) = 2ft, MSt). 

Let {D„...,D m } and {Di,... t Dfi,...,Di f ,.. ^JDlr) be 
optimal schedules of the dynamic traffic for given 
[S u ...£ m ] and {S!,..-,^*,...,^,...,^}, respect- 
ively. Then 



*=1 /=! 



Proof. From Lemma 1, there exists a decomposition 
of £>*, {At,... such that D k = 2%D 1 ^ is 
a switching matrix and L(S? k ) = Z>(5* + £>i) for 
/=!,.. The schedules_{5j,... f S?i,-.- 1 Si,---^Sr} 
and {£){,..- ^>f»,..^L- therefore satisfy 
Equations (l)-(5). Since for given 

is an optimal schedule among all possible schedules 
for dynamic traffic, satisfying Equations (l)-(5), 



m 



m Ik m <*k 

ISM * ES^I = Ea 



This theorem implies that if we apply a partial 
optimization algorithm after decomposing the 
switching matrix of the static traffic, we can obtain a 
better solution for the dynamic traffic. For example, 
Figure 1 is the example presented in Reference 10. 
By applying the local-optimal algorithm to this 
example, we obtain a schedule where the total num- 
ber of packets assigned is 21. But, this number can 
be increased by first decomposing the traffic matrix 
S into five switching matrices S U ...,S S , as in 
Figure 2. The algorithm in Section 4 contains the 
method for decomposition. Notice that the length of 
the schedule is unchanged. Figure 2 also shows the 
result obtained by applying the local-optimal algor- 
ithm starting with the TDMA frame {S,,..-,S 5 }. 
Notice that the total number of packets assigned is 
increased from 21 to 33. The upper bound derived 
in the previous section for this example is 33. Hence 
the schedule in Figure 2 is actually a true optimal 
solution. In this case, decomposition and a partial 
optimization approach generates a fully optimal sol- 
ution. However, in general this is not the case. 

4. FULL OPTIMIZATION ALGORITHM 

In this section we shall propose a full optimization 
algorithm which produces the static and dynamic 
switching matrices interactively. This algorithm gen- 
erates switching matrices for the static traffic while 
considering the possibility of inserting additional 
dynamic traffic to the generated static switching 
matrix. 

We first construct a 2Nx2N static traffic matrix 
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Figure I. An example of matrix S and matrix D 
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Figure 2. A schedule for the example in Figure 1 : the entries of D are encircled. The amount of dynamic traffic scheduled is 33 



5 from S y such that the total traffic in each line of 
3> is equal to !T max , by augmenting the dummy traffic 
matrices X 9 Y,Z to S as follows: 



(S 
Y 



3- 



where X,Y are diagonal matrices. Elements in X,Y 
are determined by the condition that the total traffic 
in each line of S is equal to r max . Let Z be a non- 
negative jnatrix such that the total traffic in each 
line of S is equal to T mBx . By augmenting the 
dummy traffic matrices, the total static traffic in 
each line of S is equal to T max , resulting in a total 
static traffic to be scheduled 2NT mj ^. It is known 
that an optimal TDMA frame for S with a length 
equal to T max can be easily generated. 9 

The matrix S can be represented as a bipartite 
graph G^UyVJE), where the nodes in U correspond 
to the rows of S, and the nodes in V to the columns. 
The edge set E contains an edge joining the nodes 
u, and v p denoted (w,,v y ), if and only if, w 0 0. A 
numerical value, ~5 ip is assigned to each edge. The 
resulting graph will be called the network associated 
with S. For this graph, we find a complete matching. 
The set of cells in S corresponding to this complete 
matching, provide a switching matrix S x for the 
static traffic. Suppose that a matched edge (u h Vj) 
corresponds to dummy traffic. This dummy traffic 
can belong to one of three matrices X,Y,Z. If it 
belongs to X, then 1 < i < N, j = N + /', and row i 
in S x has no non-zero elements, i.e. row / is exposed. 
The dynamic traffic can be assigned to row /' and 
"Sij^i is an upper bound on the maximum schedulable 
dynamic traffic for row /. Similarly, if it belongs to 
then 1 <y < N, i = TV +y, and column j is 
exposed. The dynamic traffic can be assigned to 
column j and ~S Nyj is an upper bound on the 
maximum schedulable dynamic traffic for column j. 
If it belongs to Z, then W+ 1 < /, j < 27V. Since we 
are interested in the assignment to rows 1,...^V and 
columns we simply neglect this case. 



To assign the dynamic traffic, we perform an 
assignment procedure. Let DR,(DCj) be the sum of 
row 1 (column j) of matrix D and w£ = y iVV+J . — 
DR t {wf = 3„ +JJ -DCj) for ij= 1,2,...^. If m*(m£) > 
0, then vi^(m^) represents the minimum possible 
number of unused surplus slots in row i (column j) 
even after the assignment of dynamic traffic. If 

rK,-y rvy y ^ liiwxi ou biupiub biuU) tuuiu l/v uovu anu 

the minimum possible number of unused surplus 
slots is zero. 

The dynamic traffic with the smallest row and 
column surplus slots may cause a bottleneck in the 
scheduling. Hence w 0 r = max{wf,0} + max{w?,0} is 
used as the weight for non-zero entry d tJ for 
ij = 1,2,...^V. Since the dynamic traffic can be 
assigned to the exposed lines, the calculation of w 0 
is executed only for these exposed lines. Therefore, 
if both row /" and column J are exposed, we calculate 
Wg. Otherwise, we set a sufficiently large number 
to Wij, and then solve the weighted matching prob- 
lem for the dynamic traffic which minimizes 
XiltjWjj. The cells in D corresponding to a solution 
of this weighted matching problem are candidates 
for the switching matrix for dynamic traffic. Among 
these cells, only the cells corresponding to the 
exposed rows of Si can be scheduled, since a switch- 
ing matrix can have at most one non-zero element 
in each line. Let D y be a switching matrix for the 
dynamic traffic constructed in this way. 

The matrix S x and D x are therefore the first set 
of switching matrices for the static and dynamic 
traffic, respectively. However, it would be more 
efficient to truncate these matrices so that all non- 
zero entries are equal. Let 5, be the smallest non- 
zero value of cells in ~S, corresponding to the set 
of matched edges. Similarly, let 5 2 be the smallest 
non-zero value of cells in D, corresponding to the 
set of weighted matched edges. Let 8 = min{5,,8 2 }. 
Then truncate all non-zero elements in S x and £>, 
to the value 5. This completes the first iteration for 
the generation of S x and Z),. Subtract S x and £>, 
from S and £>, respectively, and then repeat. We 
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can easily show that the switching matrices con- 
structed in this way satisfy Equations (l)-( 5). 

Algorithm Time slot assignment for heterogeneous 
environment (TSAH) 

Step 0: Initialization 

Let S={Sij} and D={d 0 } be given traffic mat- 
rices; set k 1 ; 

Step 1: 

1.0. (r^j 

SRj = 2£,j iy , i = 1,...^V; 

r max = max lV {5^C jr }; 

1.1. Create augmented traffic matrix 
S«-0; 

y,y < ^ y, ij — 1 , - • • ,7V , 

* ^max ~ SCj,j = l,...,Af; 

For i = l,...,Af and y" = 1,...^V, 

if SRtf+i (= S^iS^H-,^) — SC/y+y 

f — X2Nyr \ *U„„ -m HP CD 

V - ^i=W»VV+/V ""-xi *hh-ijf+j *^ a max — 

else J^/+ivv+y * T max — SC^ + y; 

Ste/? 2: (Form switching matrices) 

2.0. — 0; Z)* — 0; 

2.1. Set up the network associated with 5; 

2.2. Find a complete matching M in the network 
associated with 5; 

Af = {(u u m(u x )\(u 2 ,m(u 2 )\ ...^^M^))}, 
where /»(«,-) is the mate of w,; 
R — 0; C — 0; 

For w, = 1,...^V, if > then R <— /? 

+ {",}; 

For u { = N + l,-..,2Af, if m(w,) < TV then C 

2.3. Create an iVxN weight matrix W= {w 0 } 
as follows: 

For i e /?, if DR f (= SJH,^) > 0 then w; 
= J lVV+/ - DRs if wf < 0 then v< <— 0; For 
y e C, if DC/= 2£,</y) > 0 then = 
3 N+JJ - DC/, if vv? < 0 then Wj «— 0; For ij 
= \, - yN, if i e e C, and d u ^ 0 then 
w tf = + u^, otherwise, w, y = oo; Find a 
weighted matching A which minimizes 
2,2,*',, in the network associated with W\ 
A = {{x u a{x x )Ux 2 ,a{x2)\'lx N ,a(x N ))}\ 
where <?(*,) is the mate of jc f ; 

2.4. For i = 1,. . .^V, if i g R then store s J>I(/) in the 
matrix S k , else store d iM0 in the matrix 

2.5 8, = min{y Mltm(l j,...^^^}; 

o 2 = mm{d x . Mx .y\d XiXiiXi) > 0, i e { 1,. ..,//} }; 
8 = min{5,,5 2 }; 

Form a switching matrix + by replacing 
all its nonzero entries by the value 6; 



Step 3: (Updating and termination) 
S <— S - S ky D <— D - D k \ 
If S contains no non-zero entry, then STOP. 
Otherwise, set A' — k + 1; go to Step 1. 

The computation of a complete matching in 
Step 2.2. and the weighted matching computation in 
Step 2.3. can be carried out in 0(N*) time. 11 Since 
Step 2 is executed at most r max times, the overall 
worst-case time complexity is 0(r max A^). 

5. SIMULATION RESULTS 

TSAH is implemented in Pascal and simulation tests 
have been performed. Traffic matrices of 4 x 4 and 
7x7 are tested. For each size, dense and sparse 
matrix cases are tested. Dense cases have few null 
entries and sparse cases have nearly 50 per cent null 
entries. The static traffic matrix entries are randomly 
generated from the uniform distribution between 
[QyP], where P is the largest possible entry in S. 
The dynamic traffic matrix entries are randomly 
generated from a Poisson distribution with a mean 
of 4-5. For each case of matrix size and density, 
four different values of P are tested. For each test, 
we generate 200 problems and report the average. 

Tables I— IV show the results of our simulation 
experiment. They present the upper bound obtained 
in Section 2, the number of dynamic traffic sched- 
uled and the average percentage deviation from the 
upper bound for each algorithm. In Tables I-IV, LO- 
GG is the local-optimal algorithm in Reference 10, 
combined with the algorithm in Reference 8, which 
generates a TDMA frame for the static traffic. The 
algorithm in Reference 8 generates switching matr- 
ices with equal non-zero entries, which are a quite 
decomposed TDMA frame. As we have shown in 
Section 3, the local-optimal algorithm generates a 
better solution for a decomposed case. If the local- 
optimal algorithm is combined with other algorithms 
which generate a TDMA frame with less decom- 
posed switching matrices, the results of the local- 
optimal algorithm would be much worse. 

Tables I-IV show that the solution quality of our 
algorithm TSAH, is much better than that of the 
local-optimal algorithm, even though the local-opti- 
mal algorithm is combined with a decomposing 
algorithm in Reference 8. Notice that the percentage 
deviation is measured from the upper bound, UB, 
not from the true optimal value. We think the results 
of TSAH are very close to optimal values. 

As the value of P increases, i.e. as the difference 
between the quantity of static traffic and the quantity 
of dynamic traffic gets larger, TSAH provides better 
solutions. Tables I-IV also show that TSAH has 
good performance for sparse matrices. 

6. CONCLUSIONS 

In this paper, we considered the time slot assignment 
for a heterogeneous environment in which circuit- 
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Table I. Number of dynamic traffic scheduled 



4x4 Dense Matrices: 



Algorithm 


P = 10 


P = 20 


P = 40 


P= 100 


UB 

LO-GG 
TSAH 


28 

26 04 (9-43) 
26-65 (7-51) 


37 

32-24 (14-66) 
34- 14 ( 9-93) 


44 

37-71 (15-86) 
4101 ( 8-58) 


47 

39-32 (17 07) 
43-38 ( 8-76) 


Numbers in 


parentheses are the average percentage deviation from the upper bound. 




Table II. Number of dynamic traffic scheduled 




4x4 Sparse Matrices: 


Algorithm 


P = 10 


P=20 


P = 40 


P= 100 


UB 

LO-GG 
TSAH 


17 

15- 94 (11-59) 

16- 65 ( 7-97) 


22 

19-98 (10-28) 
21-20 ( 5-44) 


24 

21-84 ( 9-73) 
23-48 ( 3-61) 


23 

20-71 (12-20) 
22-58 ( 4-93) 



Numbers in parentheses are the average percentage deviation from the upper bound. 



Table III. Number of dynamic traffic scheduled 



7x7 Dense Matrices: 


Algorithm 


P = 10 


P = 20 


P = 40 


P= 100 


UB 


92 


132 


154 


167 


LO-GG 


85-53 (8*49) 


116-41 (12-65) 


131-55 (15-23) 


142-63 (15-21) 


TSAH 


86-09 (7-90) 


119-71 (10-35) 


138-80 (10-63) 


151-26 (10 07) 



Numbers in parentheses are the average percentage deviation from the upper bound. 



Table IV. Number of dynamic traffic scheduled 



7x7 Sparse Matrices: 


Algorithm 


P= 10 


P = 20 


P = 40 


/>= 100 


UB 


64 


76 


84 


85 


LO-GG 


59-30 (9 05) 


70-47 (7-91) 


77-95 (7-66) 


78-52 (8-42) 


TSAH 


60-41 (7-38) 


72-65 (5 15) 


81-33 (3-72) 


82-48 (4-00) 



Numbers in parentheses are the average percentage deviation from the upper bound. 



and packet-switched traffic share the same satellite. 
Our algorithm provides solutions better than the 
existing algorithm and has relatively low compu- 
tational complexity. 

Our algorithm is also applicable to a case where 
the traffic is homogeneous in nature and has differ- 
ent priorities. A reasonable approach may be to 
schedule the high priority traffic as the static traffic 
and the low priority traffic as the dynamic traffic. 
Extending our algorithm to a case where the traffic 
has more than two different priority types would be 
an interesting study for future research work. 
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